%T2
%二维插值
clc;clear;
N=10;%插值个数

h=[6.36,6.97,6.23,4.77;
    6.98,7.12,6.31,4.78;
    6.83,6.73,5.99,4.12;
    6.61,6.25,5.53,3.34];
x=1:4;
y=1:4;
%[x,y]=meshgrid(1:1:4);%产生坐标数据矩阵

%% 网格节点插值
method=char('nearest','linear','spline','cubic');%插值方法
x1=1:(4-1)/N:4;
y1=1:(4-1)/N:4;

s(1)=subplot(321),mesh(x,y,h),title('原网格图')

for k=3:6
   s(k)=subplot(3,2,k);%图形分割
   h1=interp2(x',y,h,x1',y1,method(k-2,:));%二维插值  
   % !!!注意!!!这里必须要把x和x1转置一下才能获得矩阵h，否则得到的h行列数不对!
   mesh(x1,y1,h1);title([method(k-2,:),'型插值曲面'])
end

axis([1 4 1 4 0 10])%设置坐标轴范围

pause
%% 散乱节点插值
x2=1:(4-1)/N:4;
y2=1:(4-1)/N:4;

subplot(2,1,1);
mesh(x,y,h),title('原网格图')

h2=griddata(x',y,h,x2',y2);%随机点插值
subplot(2,1,2);
mesh(x2,y2,h2);
title('散乱数据点插值曲面')







